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We introduce the interaction cost of a non-local gate as 
the minimal time of interaction required to perform the gate 
when assisting the process with fast local unitàries. This cost, 
of interest both in the areas of quantum control and quantum 
information, depends on the specific interaction, and allows to 
compare in an operationally meaningful manner any two non- 
local gates. In the case of a two-qubit system, an analytical 
expression for the interaction cost of any unitary operation 
given any coupling Hamiltonian is obtained. One gate may 
be more time-consuming than another for any possible inter- 
action. This defines a partial order structure in the set of 
non-local gates, that compares their degree of non-locality. 
We analytically characterize this partial order in a region of 
the set of two-qubit gates. 
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An elementary concern in quantum information the- 
ory is to stablish the trade-off between different physi- 
cal resources that are relevant for information process- 
ing. A controllcd Hamiltonian interaction between quan- 
tum systems is one instance of useful resource. It can be 
employed, for example, to simulate the dynamics of an- 
other multipartite quantum system. On the other hand 
multi-particle unitary gates are a requirement for univer- 
sal quantum computation. In particular, two-qubit gates 
— together with one-qubit gates — can be taken as the 
building block of quantum computers. 

A detailed study of the connections cxisting between 
non-local Hamiltonians and non-local gates is thus of in- 
terest from a quantum information perspective, but this 
issue is also relevant in other areas. For instance, the 
synthesis of multipartite gates from Hamiltonian inter- 
actions — and, in particular, time-minimizing schemes — 
has been recently analyzed in the context of quantum 
control theory [0. Whereas the requirements for arbi- 
trary manipulation of single qubits are presently met in a 
number of experimental schemes, the enginecring of two- 
qubit gates can be only (partially) achieved with very few 
systems Q. In real experiments not only an interaction 
Hamiltonian between the qubits, but also considerable 
command on them in order to process the interaction, are 
required. For instance, mechanisms to switch on and off 
the interaction, as well as to accurately drive the systems 
towards the desired joint evolution, are needed. But even 
from a theoretical perspective, a description of two-qubit 
gates in terms of interactions able to prescribe optimal 
protocols for gate synthesis was so far missing. Here we 



shall provide such a description. 

More generally, we consider a set of subsystems with 
a given Hamiltonian H, and assume that arbitrarily fast 
local unitàries (LU) can be performed to properly tailor 
the evolution that H induces. The aim is to perform some 
joint unitary transformation U on the systems. This is 
the setting considered in jïj and corresponds to the so- 
called gate simulation under LU of ||. Two definitions 
are needed to specify the problems that we shall address. 

Definition 1: The interaction cost Ch(U) of a non- 
local gate U given a Hamiltonian H denotes the minimal 
time needed in order to perform U using the interaction 
H and fast LU. 

Definition 2: We say gate U is more non-local than 
gate V, and write V < U, when for all interactions H 
the interaction cost of U is never smaller than that of V, 



V < U = C H (V) < C H {U) ViT. 



(1) 



First we shall show how the interaction cost Ch can be 
explicitly computed for any gate and any interaction of 
a two-qubit system. This is possible by considering re- 
sults recently developed in the areas of quantum control 
jïj and quantum information In considerable 

progress towards the solution was made, and only a final 
optimization was left unsolved. The results of pro- 
vide the tools needed to perform such an optimization 
and thereby complete the results of ||]. 

Definition 2 introduces a partial order structure in the 
set of non-local gates. This structure captures the intu- 
ition, in terms of the resources needed to perform a gate, 
that one gate may be "more non-local" than another. 
Our second result is an analytical characterization of this 
structure in a region of the set of two-qubit gates. 

We start by describing known facts concerning the sim- 
ulation of non-local Hamiltonians and the synthesis of 
non-local gates. 

(i) Optimal simulation of two-qubit Hamiltonians un- 
der LU. Any Hamiltonian acting on two qubits is 
uniquely represented, for the purposes of simulation un- 
der LU, by its canonical form pipí 



H =} hiOi 



h\ > h 2 > \h 3 



(2) 



where o~i, i = 1,2,3, stand for the Pauli matrices. In 
the rest of the paper H denotes a Hamiltonian written in 
its canonical form and h denotes the vector (hi, /12, 
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with its properly ordered coefficients. The special ma- 
jorization relation x -< s ybetween three dimensional real 
vectors x and y is relevant in this context. It is given by 
the set of incqualitics 

x\ + x 2 - x 3 < yi + y 2 - 2/3, 

x\ + x 2 + x 3 < 2/1 + 2/2 + 2/3- (3) 

where the components Xi and yi are assumed to fulfi.il 
xi > x 2 > \x 3 \ and y 1 >y 2 > \y s \ 0. 

Fact 1 (Theorem of Q): The minimal time over- 
head tH'\H the inverse of efhciency sh'\h) required 
to simulate Hamiltonian H' by Hamiltonian H and fast 
LU is the minimal valuc of c > such that the vectors h! 
and h satisfy h! -<. s ch. Protocols for optimal simulation 
are known. 

(m) Optimal synthesis of two-qubit gates under LU. 
Any two-qubit gate U can be written in terms of local 
unitàries Uj and Wj and Pauli matrices as 0J^] 

U ={u 1 ®v 1 )e- l ^* Xkak ® ak {u 2 ®v 2 ). (4) 

In |^| it is shown how to obtain this decomposition. No- 
tice that gate U is equivalent, up to local unitàries per- 
formed on the qubits before and after U , to 

[T^e-íE*** 0- *®**. (5) 

Since we assume the ability to perform instantaneous 
(i.e., sufhciently fast) LU operations, the synthesis of U 
is as time-consuming as that of Ur, and we need only 
focus on the later. In addition ||, to each U there cor- 
responds a unique U? with A° = (A^A^Ag) such that 
A? > A° > |A°|, A°,A° e [O.tt/4] and A° e (-*r/4,7r/4], 
that we shall call its canonical form. In what follows we 
will often represent any two-qubit gate U by its canonical 
form U Xo or by its corresponding (unique) vector A° = 
(Xi, X 2 , X 3 ). Recali that all commutators \a 3 -®a d ■, (Jk®<Jk\ 
vanish, and that cxp(— m/2o~j ® o~j) = —io-j®<jj is a local 
gate. This implies that for any vector n — (ni,n 2 ,n 3 ) 
with integer components rij, 

U Xo = LU U Xo U* Pl = U Xo+in , (6) 

with A° + 7r/2n essentially exhausting all vectors com- 
patible with the gate U x<> (|]. [—lu is used to denote 
equivalence under LU.]. 

In Theorem 10 of Q| the problem of time-optimally 
producing a two-qubit gate U using interaction H is 
shown to reduce to a specific minimization over all possi- 
ble decompositions of U of the form . Here we rephrase 
the theorem in terms of the notion of Hamiltonian simu- 
lation and the concepts introduced before. Without loss 
of generality, we refer only to unitary operations that can 
be written as in (|^), and associate a self-adjoint operator 
H x = J2i Xí<Jí <8> <Ji to each possible decomposition. 



Fact 2 (Theorem 10 of [0, readapted): The time- 
optimal way to synthesize gate U with interaction H and 
fast LU consists of simulating, among all Hamiltonians 
Ht such that U — exp(—iH x ), the one with smallest 
time overhead t H ^ H . The minimal interaction time (i.e., 
the interaction cost Ch(U)) is given by the smallest time 
overhead íh^h- 

Our first aim is to perform the optimization described 
in Fact 2. This is feasible because we have an analytical 
characterization both of all (infinitcly many) decomposi- 
tions of U (cf. Eq. (^])) and of the time overhead s Hx \ H 
for any decomposition (cf. Fact 1), as expressed in the 
following lemma. 

Lemma: The interaction cost Ch(U) is the minimal 
value c > such that a vector n of integers exists satis- 
fying 

A° + |n < s ch. (7) 

It is useful to introduce, for each n, the pre-cost Cj% as 
the minimal value c > such that A° + 7r/2n -< s ch. Pre- 
cost ca is the overhead needed to simulate H Xo+7r ^ 2fí by 
H or, equivalently, the minimal time t needed to travel, 
in the set of non-local gates, from the identity operator 
to U along the path defined by A° + -ïï/2n. Intuitively, 
a large n corresponds to a "long" — and therefore non- 
optimal — path. Following this intuition we arrive at our 
first result. 

Theorem 1: The interaction cost Ch(U) or minimal 
time needed to create gate U by using Hamiltonian H 
and fast LU is given by 

C H (U) =min{ 0(0,0,0), c (-i,o,o)}, ( 8 ) 

that is, the minimal of two pre-costs, one correspond- 
ing to the canonical vector A° = {Xf, X 2 , X 3 ) of U and 
the other to the vector (| - A?, A§, -A§) [§. The time- 
optimal protocol consists in simulating the correspond- 
ing Hamiltonian (either h\ = (A^A^A^) or h 2 = (? — 
A", A°, — Ag)) by H for time t = C H (U). 

Remark. Thus, in order to time-optimally perform 
gate U with Hamiltonian H, we can proceed as follows. 
Using Ref. ||, we compute A° from U, and using Re fs. 

we compute h from H . Theorem 1 gives the minimal 
time of simulation and the Hamiltonian (either h\ or h 2 ) 
to be simulated, and finally Ref. ||] describes an optimal 
protocol for simulating the convenient Hamiltonian by H 
and LU. 

Proof: We need to see that Ch{U) as given by Eq. 
(||) is the minimal pre-cost, i.e. Ch(U) < Cfj for all 
n. It is straightforward to check from Eq. (||) that 
(Tl.i) for any two vectors x and y, with components 
x\ > x 2 > \xs\, 2/i > 2/2 > I2/3I, the minimal c > such 
that x ^ s cy satisfies c < 3xi/yi; (Tl.ü) if x ^ s x\ 
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then x' -< s y x -< s y {-< s is a partial order!). In 
particular, let d > be the minimal value such that 
x' -< s dy. Then x < s x! =}■ x -< s dy, so that the min- 
imal c > such that x -< s cy always satisfies c < d. 
Now, recali that by definition tt/4 > A° > 0, and no- 
tice that if some component nj of n fulfills \nj\ > 1, 
then the maximal component of the reordered version 
of A° + 7r/2n is at least 3ir/4. Thus, because of (Tl.i), 
A° ~< s A° +ir/2ïí. Then (TI.m) implies that c^o.o) < c fi- 
Therefore we can restrict our attention to vectors n with 
\rij\ < 1. A case by case check shows that the pre-costs 
with n e {(-1,-1, -1), (0, -1, 0), (0, 0, -1), (0, 0, 1)} 
fulfill C(_ 10 ) < Cfi, since (cf. point (TI.m) above) 

A°+7r/2(-l,0,0) -< s X q + tt/2tÍ fo|. Similarly, we obtain 
that for the remaining vectors n with \rij\ < 1 the pre- 
costs satisfy C(o,o,o) < ca, because A° + 7r/2(0, 0, 0) -< s 
A° + ir/2n. The only remaining configurations, with vec- 
tors n 6 {(— 1, 0, 0), (0, 0, 0)}, are incomparable accord- 
ing to the -< s relation — unlcss A? + |A§| < 7r/4, in which 
case we always obtain C(o,o,o) < c (-i.o,o) — i an d this is 
why the optimization of Eq. (^|) has to be performed. □ 

Corollary: (a) When U is such that A? + |A° < tt/4, 
then the interaction cost is always given by 

C H (U) = c^o.o)- (9) 
(ò) If, instead, Aj" + 1 A§ | < 7r/4, then Hamiltonians H and 
H' always exist such that Ch{U) — C(o o o) < c (-ioo) 



Ji' 



h' 



and C W (U) =c' ( L 1A0) < c (0 0) . 

Proof: (a) follows from the fact that A° + | A§ | < 
tt/4 => A° -< s A° + tt/2(-1,0,0). (b) can be checked 
by considering H and H' given by h = A° and h! — 
(7r/2-A?,A^-A°). □ 

In order to analyze Eq. (||) we first consider some ex- 
amples. For the Ising interaction H = ha^ ® 03 (equiva- 
lently ho\ <g> <j\) and an arbitrary gate U, Eq. (g) reads 
(cf. Theorem 2 of @), 



A? + A° + | A§ | 



(10) 



Let us now instead focus on three specific gates and arbi- 
trary interactions. By |m) (gi |n) (m, n = 0, 1) we denote 
the computational basis of two-qubits. The CNOT gate 
is dcfincd as 



|m) ® |n) — > \m) ® \n © m), 



where © is sum modulo 2. Using the method described 
in Ref. [|| we obtain its canonical vector, A° — {tt/4:, 0, 0). 
Similarly, the SWAP gate, 



|m) <g) |n) 



|n) (g) |m) 



(12) 



has vector A° = (n/4, tt/4, tt/4). We also consider a third, 
intermediate gate Uxy with A° = (tt/4, tt/4, 0), that cor- 
responds to 



\m) ® \n) — » il m -"l|n) ® |m). 
For these three gates we find 



Cíï(CNOT) 
Ch{Uxy) 
C H (SWAP) 



7T 1 

~4h~i 

TV 2 



4 h\ + h% - \h 3 \ 

TT 3 

Ah 1 + h 2 + \h 3 \ 



(13) 

(14) 
(15) 
(16) 



With these examples at hand we make the following two 
observations. First, to any gate U there corresponds a 
natural interaction Hjj, with vector either h\ or h 2 as 
defined in Theorem 1. This natural interaction allows 
to perform gate U optimally without need to interme- 
diately simulate another Hamiltonian and therefore the 
time inefficiencies inherent in the process of simulation 
are avoided. In this sense the natural interactions for 
the CNOT gate, gate U X y and the SWAP gate are, re- 
spectively, the Ising interaction <j\ (g> 0\ , the XY-model 
interaction o~± ® o~\ + o~2 ® &2 and the Heisenberg or ex- 
change interaction o\ ® a\ + 02 ® 02 + 03 8> C3. 

The second observation is that for any fixed interaction 
H, e.g. H = hai ®ai as in Eq. (|Ï0|), the interaction cost 
induces an order in the set of gates. For instance, accord- 
ing to Eq. (|Ï0|), a SWAP is the most time-consu ming 
gate when the Ising interaction is available. Eqs. ( [14] - 
p^ ) also show, however, that such an order depends on 
the available interaction. Using the exchange interaction 
H = a\ ® <ti + (T2 8) 02 + 03 ® (T3 , Uxy is twice as time- 
consuming as a SWAP gate. 

Let us move to Definition 2. It endows the set of non- 
local gates with a partial order structure based on the 
notion of interaction cost, but which is independent of 
any particular interaction. By comparing the resources 
required to perform two gates, such a partial order cap- 
tures the intuition that some gates are more non-local 
than othcrs. 

We have already argued that no gate more non-local 
(i.e., more time-consuming for all interactions) than all 
the others exists. It is also easy to see that a gate a\° is 
always less non-local than A° for any a € [0, 1], since the 
pre-costs are linear in a. Next we present an analytical 
characterization of the partial order relation V < U in a 
region of the set of two-qubit gates jïï| . 

Theorem 2: Let U and V be two two-qubit gates 
(11) with corresponding ordered vectors \\ and Ay such that 



in both cases the restriction AÇ H~ | A§ | < tt/4 holds. Then 
gate U is more non-local than gate V if and only if Ay -< 3 
\° 



V<U O 



-<.. 



Ajj. 



(17) 



Proof: Recali that the restrictions on A^ and Ay im- 
ply, because of Corollary 1, that the interaction costs 
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Ch(U) and Ch(V) are given, respectively, by the small- 
est Cu, cy > such that 



X v -< a c v h, 
X v -< s c v h. 



(18) 
(19) 



Suppose first that V < U , that is, that for any Hamil- 
tonian H we have Cu(V) < Ch(U). Then we also have 
Xy -< s Cl·l(U)h. In particular, if we choose the inter- 
action H to have vector h = A^, we have Ch(V) < 
C H (U) = 1 and X v -< s C H (U)h = A^, which proves the 
direct implication. The inverse implication follows from 



(Tl.n) of the proof of Theorem 1, which shows X v -<, 
implies that C H (V) < C H (U) for all H. □ 

As an example of this result, we see that the Uxy 
gate is more non-local than the CNOT gate, and that, 
as it was to be expected, gates with sufhciently small 
components |A° are less non-local than those with large 
|A?|& 

Finally, recali that entanglement can be used as a cat- 
alyzer for Hamiltonian simulation ||Ï3| . In particular, if to 
each of the two interacting qubits A and B we attach an 
extra qubit A 1 and B' , where the pair A'B' is in a maxi- 
mally entangled state, Hamiltonian H between A and B 
can be used to perform, without consuming the entan- 
glement of A'B' , more powerful simulations than before, 
provided that fast LU are allowed in AA' and BB' . Con- 
sequently, the interaction cost of gates is modificd when 
not only LU, but also entangled ancillas are available. 

In this work we have characterized the time-optimal 
synthesis of two-qubit unitary transformations using an 
arbitrary two-qubit Hamiltonian. In particular, the in- 
teraction cost Ch(H) has been computed and optimal 
protocols have been described. We have also character- 
ized, in a region of the space of two-qubit gates, a partial 
order structure related to their degree of non-locality. 
These results can be applied to the study of the interac- 
tion cost for particular processes, such as the creation of 
a maximally entangled state or the transmition of 
a classical or quantum bit of information from one qubit 
to another Q. All these discussions involve only two 
interacting qubits. It would be desirable to obtain a gen- 
eralization to highcr-dimcnsional systems. The lack of 
an analog to decomposition (|4|) in these cases is a serious 
drawback. Another interesting generalization consists in 
considering the asymptotic scenario, where the aim is to 
perform a large number of copies of the same gate. 
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The s-majorization relation can be extended to non- 
ordered vectors as follows. Given a vector ü = 
(ui,U2,ug), we construct a new vector Ü" = (ui, u^u^), 
tíf > tí| > |tí.§ I by first decreasingly reordering the modu- 
lus of the components Ui, and by then giving 113 the sign 
of the product U1U2U3. Then, for any pair of vectors ü 
and v, u ~< s v denotes the set of inequalities (^) applied 
to u 3 and v" . 

[8] By permuting the components of A° + Tv/2n, as well as 
multiplying by — 1 two of the components, other vectors 
can be build that also correspond to a decomposition 
of U. These are, nevertheless, equivalent to the original 
decomposition from the perspective of the s-majorization 
relation -< s , and need not be analyzed for the purposes 
of optimal gate synthesis. 

[9] (f - Aj.A^-As) is the reordered version @ of A° + 
7r/2(-l,0,0). 

[10] Let us consider, as an example, n = (0, 0, 1). If we prop- 
erly reorder j?j the coefRcients, we obtain the vector 
(7r/2-A ( í,A^,-Ag) from A°+ tt/2(- 1,0,0) and the vector 
(tt/2 + X° 3 , A?, X° 2 ) from A° + tt/2(0, 0, 1). Then 

7T/2-A? < tt/2 + A° 
71-/2 - A? + AÍ^ + A° < tt/2 + A" + A? — AÍj 
tt/2 - A? + A° - A° < tt/2 + A° + A? + A°, 

so that indeed A° + tt/2(-1, 0, 0) ^ s A° + tt/2(0,0, 1) and 
therefore the pre-costs satisfy C(_ 10 ,o) < c (o,o,i) (cf. point 
(Tl.n) in the proof of Theorem 1). 

[11] Contrary to the problem of Hamiltonian simulation, 
where the analogous partial order has been solved com- 
pletely 0], here we do not have a linear space of Hermi- 
tian operators (Hamiltonians), but the rather involved 
manifold S , Í7(4). In particular, in the region surrounding 
A° = (7r/4, 7r/4, 7t/4) geometrical aspects of this manifold 
make the characterization very difHcult. We have already 
met a manifestation of this fact in Eq. ([J), where two pre- 
costs, instead of just one, have to be considered due to 
periodicity properties of unitary evolutions. The analyt- 
ical characterization presented here is possible precisely 
away from A° = (7r/4, 7t/4, 7t/4). 

[12] Remarkably enough, the necessary and sufficient condi- 
tions for this partial order are completely analogous to 
those for emcient Hamiltonian simulation Q . Such a cor- 
respondence was, in principle, only to be expected for 
infinitesimal gates. It also holds for finite gates due to 
some convexity properties explored in Theorems 9 and 
10of§. ' 

[13] G. Vidal and J.I. Cirac, quant-ph/0108077 . 
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